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We study the effects of tfiermal fluctuations of gluons and tlie diquark pairing field on the 
superconducting-to-normal state phase transition in a three-flavor color superconductor, using the 
Ginzburg-Landau free energy. At high baryon densities, where the system is a type 1 supercon- 
ductor, gluonic fluctuations, which dominate over diquark fluctuations, induce a cubic term in the 
. - . Ginzburg-Landau free energy, as well as large corrections to quadratic and quartic terms of the 

\l ' order parameter. The cubic term leads to a relatively strong first order transition, in contrast with 

I the very weak first order transitions in metallic type 1 superconductors. The strength of the first 

. order transition decreases with increasing baryon density. In addition gluonic fiuctuations lower 

the critical temperature of the first order transition. We derive explicit formulas for the critical 
I^J^' temperature and the discontinuity of the order parameter at the critical point. The validity of the 

'~{ , first order transition obtained in the one-loop approximation is also examined by estimating the size 

of the critical region. 
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I. INTRODUCTION 
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Tjj- . Degenerate quark matter at high baryon density is expected to undergo a phase transition to a color superconducting 
■ state 0. The properties of color superconductors have been studied so far in the weak coupling regime with one- 
OA , gluon exchange between quarks [iJMi in the strong coupling regime with an effective four-fermion interaction M, 
' ' in Ginzburg-Landau (GL) theory and from the perspective of sum rules and phenomenological equations Q- 

[ A major difference of color superconductors and metallic superconductors is that the former is a Inghly relativistic 
■"".^ ■ system in which the long-range magnetic interaction (dynamically screened only by Landau damping [g ) is responsible 
' for the formation of the superconducting gap. The dynamically screened interaction leads to a nonstandard form of 
the gap, A{g) cx fj,exp{—3'!T^ / V2g) in the weak couphng region, with g the strong couphng constant and fj, the baryon 
O i' chemical potential Q • Despite this non-BCS feature of color superconductivity, the transition to the normal phase at 
D finite temperature T in mean-field theory is second order, with a BCS critical temperature Tc ~ 0.57A 'lO]. 

In this paper we address the question of the modification of the structure of the transition due to inclusion, beyond 
mean- field theory, of thermal fluctuations of the gluons and of the diquark pairing condensate ~ {tpip)- The effects 
of thermal fluctuations on the phase transition were first studied in BCS superconductors in metals in Ref. Ill, and 
. , in finite temperature field theory in Ref. [T^ . The GL free energy for a metallic superconductor without photon 
• degrees of freedom has a global U{1) symmetry. It therefore falls into the 0(2) universality class and shows a second 
" ■ ' order phase transition [ij]. However, the coupling to photon fluctuations may lead to a flrst order phase transition 
[nl; in particular, type I materials have a ^Icafc first order transition, characterized by a cubic term of the order 
parameter in the GL free energy. In color superconductors, a first order transition is also expected 2]. However, there 
are crucial differences from metallic superconductors. Firstly, the fluctuations of the diquark field alone may lead to 
a first order phase transition; in fact, the GL free energy in color superconductivity (without gluons) has a global 
color-flavor symmetry, SUc{3) x S'C/l+j?X3), which exhibits no infrared flxed-point in the renormalization group flow 
of the coupling constants, and is thus likely to show a first order phase transition 1^ Il5j|. Secondly, thermal gluon 
fluctuations may induce a relatively strong first order transition, in contrast to the metallic case, partly because of 
the relativistic nature of the quarks and partly because of the large coupling constant as = g 

We study the effects of fluctuations of the pairing and gluon flelds on the phase transition via their effects on the GL 
free energy, emphasizing the relative importance of the diquark and gluon fluctuations, and the need to treat gluon 
fluctuations consistently, keeping all terms of the same order. We estimate, semi-quantitatively, the strength of first 
order transition as well as the modification of the transition temperature. Unlike the conclusion of for two-flavor 



The color superconducting tran sition is discussed from the point of view of the Thouless argument on fluctuations of the pairing field 
in the normal state in Ref. l25l . 
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color superconductivity, we find that the first order transition weakens with increasing baryon density, and that the 
transition temperature is lowered from its mean-field value. 

In Sec. II, we review the GL approach to color superconductivity, following 0,01 ■ We consider two characteristic 
pairings in three-flavor superconductors: color-flavor locking (CFL) and isoscalar (IS) ordering. Then we discuss 
the size of the thermal fluctuations of the pairing field and the gluons. We study the relative magnitude of these 
fiuctuations and the question of whether the system is type I or type II in the framework of GL theory for a three-flavor 
color superconductor. The validity of the one-loop approximation to evaluate gluon fluctuations is discussed for the 
type I case. In Sec. Ill, we focus on type I superconductors realized in the weak coupling regime, and calculate the 
effects of the thermal fluctuations of the gluons in the one-loop approximation. A first order transition is induced for 
both CFL and IS pairings. The strength of the transition and the critical temperature are evaluated explicitly for the 
CFL state at high density. The proposed critical end point of the first order transition in the low density regime 
is beyond the scope of this paper and will not be discussed here. Section IV is devoted to a summary and concluding 
remarks. In the Appendix, we summarize the parameters of the GL free energy in the weak coupling regime Q. 



II. GINZBURG-LANDAU FREE ENERGY 



As a prelude to analyzing the color superconducting phase transition at finite temperature, we first review the GL 
free energy and the pairing fields in the absence of fiuctuations. We then estimate the size of fluctuations around the 
mean-field and the critical regions both for gluon and diquark fields within the Gaussian approximation. 



A. Three-dimensional effective theory 



Let us consider a system of degenerate massless u, d, s quarks with a common Fermi momentum. The pairing 
gap of a quark of color b and flavor j with that of color c and flavor k in the = 0'*' channel is written as 
4'bcjk] further assuming that the pairing takes place in the color-flavor antisymmetric channel, which is expected 
to be the most attractive in the weak coupling, the gap is parametrized as 4>hcjk — ^abc^ijk da 00. Under G = 
SU{3)c X SU{3)l+r X U{1)b, the order parameter transforms as a vector, and by construction belongs to the 
(3*, 3*) representation of SU{3)c and SU{3)l+r- We consider only Cooper pairing of even parity in the present 
analysis. This is because the presence of instantons favors even rather than odd parity pairing In the absence of 
instantons, the state of even parity would be degenerate with that of odd parity, giving rise to an extended form of 
the order parameter [T^ . 

The GL free energy density in three spatial dimensions, written in terms of the order parameter field d^(x), with 
coupling to the SU{3)c gluon gauge fields, reads 

S = aJ2\da\' + f3i{J2\d,\'y + /32 E ' d.P + 2^^^ l(Ad)aP + ^GLGL. (1) 

a a ah a 

The parameters a, /3i, and (^2 characterize the homogeneous part of the free energy, while kt is the stiffness parameter, 
controlling spatial variations of the order parameter. Since is antisymmetric in color space, the color-covariant 
derivative Di is 

{Did)a=dlda+'-gAr{X^*d)a, (2) 

where the A"* are the complex conjugates of the Gell-Mann matrices, and Gim is the spatial part of the gluon 
field-strength tensor, 

G^„ = diAZ, - d^Af + gf^p^AfAZ^. (3) 

The free energy density Eq. (Q) may be interpreted as an SU{3)c x SU{3)f scalar field theory coupled to an SU{3)c 
gauge field in three spatial dimensions. Equation is model independent and valid near the critical temperature of 
the second order transition when the average value of d^ (x) is small. 

Although the general analysis does not require specific values of the parameters in the GL free energy, it is useful 
to bear in mind their characteristic scales, as found in weak coupling (see the Appendix), a ~ ln(T/Tc) and 
01 ^ P2 = 3kt ^ {p/TcY , with /J, the baryon chemical potential and Tc the weak coupling critical temperature. 

The order parameters for color-fiavor locking (CFL) and isoscalar (IS) ordering in three-flavor matter are 

d Sa^ (CFL) 

dda3S,3 (IS). 
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We shall in general take d to be real, and set the direction 3 to be & in color and s in flavor. If we consider only the 
uniform field configurations that minimize S, Eq. and neglect thermal fluctuations around the stationary value 
(the mean- field approximation), we obtain A" = A = and 

for (32 >0 and 3/3i + (32 > 

(5) 

for -(3i <(32<0. 

Since a changes sign at the mean field Tc, it is useful to rewrite it in the form, 

a = aot, (6) 

where t — (T — Tc)/Tc is the reduced temperature, and ao > 0. As is evident from Eq. (jsj, the system undergoes a 
second order phase transition from the paired state to a normal state quark-gluon plasma as T increases. Whether 
the paired state just below Tc is CFL or IS depends on the values of (3i and ^2- In the weak coupling limit where 
(3i — P2 and ag ^ /i^, CFL ordering is favored, with dcFL ^ Tc^/\t\ for T ^ Tc- 



"CFL 



a ( /5i + ift 

(32)-' 



B. Fluctuations about mean field and the critical region 



Let us now consider the effect of thermal fluctuations of the smn-zero diquark (scalar) field and the spin-one gluon 
fields about their mean values in the Gaussian approximation [ig. In this approximation, only the quadratic part of 
the fluctuations of the fields about their means (denoted by "cl"), (A")ci — and (d^)ci in Eqs. are kept in the 
free energy. The fluctuation part of the free energy is then AS = [S{d, A) — S'(dci, Aci)]quadratic- The fluctuations 
of the gauge fields at the same spatial coordinate are given by the thermal average, (A"A'^), of the product of the 
gauge fields, where 



_ /d{MKA}...e-^^»^}-{^'^>l/^ 
^'"'^ Jd{Sd}d{A}e-^Sl{A},{Sd}]/T ' 

with {A} = {A\ . . . , A^}, and {Sd} = {d-d^i} = {Re(d^) - Re(d^,i), Md^) - Md^i), 
After diagonalization of in color, we find the gluon field fluctuations, 



(A" A'') 



25a[3 T 



d^k 



k|<A 



(2^)3 k2 + (mA) 
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Ota 



(7) 

,Im(d|j)-Im(d|,,i)}- 
(8) 



where we have taken the Coulomb gauge V • A = 0. The momentum A is an ultraviolet cutoff, which corresponds 
to an upper bound on the wave numbers of the classical thermal fluctuations with zero Matsubara frequency. This 
cutoff is inversely proportional to the size of the quark pairs (A ~ d ~ Tc) In the following we take A — Tc 

for simplicity. In Eq. ©, {mA)af3 is the Meissner mass matrix, calculated in Ref. for IS and CFL orderings; the 
Meissner masses are the inverse correlation lengths of the gluon field fluctuations. The components of this matrix are 
given in Table |l| 

In weak coupling, where the system is color-flavor locked, rriA ^ 9^J'\/W\ foi' ^ ~ ^c- Since the Meissner mass 
is vanishingly small compared to A ~ Tc near the second order critical point, we can expand Eq. (jSJ in terms of 

{mA)aa/Tc: 



a \0\ — 



(A" A 



TTc I TT {■mA)aa 



ifnA)aa 

Tc 



O 



{mA)c 
Tc 



(9) 



We can similarly calculate the expectation value of the product of the fluctuations of the scalar diquark field. In 
this case it is convenient to work in the color-flavor space (a = R,G,B, i — w, d, s) in which the part of AS of 
quadratic order in 5d^ = — (d^)ci is diagonal. (For an IS condensate, the original color-flavor space provides 
the diagonalization.) For notational simplicity we write the diagonalized field as {6d)p = (JdJ^)™ where m (=1, 2) 
distinguishes the real and imaginary parts of 5d^ and the p{= a, i, m) summarizes all the indices. We then obtain 



((M)p(M).) = 



= S, 



""^^tJmkk i^^)' k^-f (m,)2 



)pp 



(10) 



Tc 



Tc 



-O 



' {md)pp ^ 
. Tc , 
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Here {md)pp is the matrix of inverse correlation lengths of the order parameter fluctuations, whose diagonal components 
are given in Table 



TABLE I: The inverse squared correlation lengths of the scalar (d) and gauge {A) fluctu- 
ations, and m\, together with the number of degenerate modes corresponding to each 
fluctuation. 
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The number of modes corresponding to a given correlation length is also indicated in Tabled The ninefold massless 
scalar modes (m^ = 0) may be understood as follows. The IS state, characterized by = d&ai^ii, Eq. is 
invariant under G' — SU{2)c x SU{2)L+ii xU{l)x U{1). Here the first C/(l) symmetry corresponds to a simultaneous 
rotation in baryon-color space, and the second to a simultaneous rotation in baryon-flavor space. Thus the number 
of Nambu-Goldstonc bosons is dim[G]— dim[G'] — 17—8 ~ 9. The CFL state, characterized by — dSai, Eq. Q), is 
invariant under G' — SU{'3)c+l+r- Thus one has 17—8 = 9 Nambu-Goldstone bosons in this case too. 

Note that not all massless scalar modes with = in Table HI are physical. Parts of them are absorbed in 
the longitudinal components of the gluon. As a result, only four massless modes out of nine are physical in the IS 
state, while only one massless mode is physical in the CFL state. For massive scalar modes, their masses behave as 
"id~Te\/^for T<Tc. 

As discussed in [llj, the initial terms in the expansions © and 1)101) proportional to TTc simply shift the critical 
temperature, Tc, of the second order transition. This is because they modify the coefficient of the quadratic term of 
the order parameter in the GL potential. On the other hand, the terms proportional to TTcirid and TTcUia in Eqs. 
and H1U|I induce a cubic term of the order parameter in the GL potential, and thus generally drive the phase transition 
to first order ■ Whether the resultant first order transition is reliable or not can be checked by estimating the size 
of the critical region on the basis of Eqs. © and T?]. The terms TTcm'^ and TTcm\ in Eqs. © and (QHI) modify 
the coefficient of quadratic term of the order parameter in the GL potential, which, as we shall see, turn out to be 
important in determining the strength of the first order transition. 

We now discuss the critical regions for scalar and gauge fluctuations. In the immediate vicinity of Tc, fluctuations of 
the soft modes become significant, leading to a breakdown of the Gaussian approximation 19]. The temperature span 
of this critical region can be determined from standard scaling arguments near the critical point (l^ . For our problem, 
the typical spatial scales of scalar and gauge field fluctuations are and m'^ , respectively. Using these scales, we 
define the "effective" coupling strengths among the soft modes for the scalar and gauge fields as fiiTd {32tt'^ n'^nid) 
and g^Tc/ {2TT^mA) — 2asTc/{nmA), respectively. We introduce the factor 27r^ since the effective coup ling s to be 
used in the perturbative expansion are always associated with the phase space factor l/(27r^) = 47r/(27r)^ These 
coupling strengths should be small enough that the calculation of the free energy in a loop expansion is meaningful. 

Also the three-dimensional effective theory for the soft modes is meaningful only when the masses of the soft modes 
are small compared to A ~ Tc- Combining the conditions discussed above, we find the necessary (but not sufficient) 
conditions for the Gaussian approximation to be valid. 



lA 



327r2K2, 



Tc < TOd < Tc. 



(11) 

(12) 



namely, that the temperature should be inside the appropriate region where the masses of the soft modes are not too 
small and not too large. Also the above equations imply that the coupling constants should be sufficiently small. 



m 



327r2K2, 



« 1, 



2a, 



< 1. 



(13) 
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FIG. 1: Schematic depiction of the relative magnitude of the critical regions, td and ea- Here Tc is the critical temperature 
for the second order transition in the mean-field approximation, and T* the critical temperature of the first order transition 
induced by gauge-field fluctuations, derived below. In the hatched region where |r — Tc] > Tc, the Ginzburg-Landau expansion 
in terms of the order parameter is not justified. 



We define the ratio of the typical spatial scales of scalar and gauge field fluctuations, 

K= ~ 2- , (14) 

which measures the relative importance of the fluctuations to the GL free energy. As we shall see in Sec. Ill, gauge 
fluctuations are more important than scalar fluctuations for k ^ 1, while scalar fluctuations are more important for 
K 3> 1. In weak coupling at high density, one finds 

oi,\ii) asVlOOMeVy V M / 

where we have neglected a numerical coefficient of order unity. Since in weak coupling, the ratio Tc//i is exponentially 
suppressed as exp(— 37r^/-\/2.g), k is considerably smaller than unity at high baryon density, and gauge fluctuations 
dominate over scalar fluctuations. 

The parameter k defined in Eq. H14I) is the Ginzburg-Landau parameter [211] that distinguishes the type of color 
superconductor under an external chromomagnetic field; k = 5 /S^ where 5 is the penetration depth, and ^ the coherence 
length. By explicit calculation of the surface energy of a domain wall separating the normal and superconducting 
phases in the presence of an external magnetic field, one finds that type I color superconductivity is realized for 
Kis < I/V2 in the IS state and for kcfl < 0.589 in the CFL state with (3i = (32 [l7j . 

As Eqs. (|14|l and H15|) indicate, rud <C ruA is realized in weak coupling. Therefore, the color superconductors are 
type I at least at very high densities. Whether they are type I or type II at low densities is not known. In the next 
section, we assume type I behavior and consider only the effect of gauge fluctuations. 

In the weak coupling limit, the Ginzburg criterion, the first of the inequalities in Eqs. Hll|l and p2|l . becomes 

, » , f^v ^ 10- r^^V f i^V. (16) 

V M / V 100 MeV / V M / 

TcV as f Tc \VlGeVy 



where we have again neglected unimportant numerical coefficients of order unity in defining ed,A- 

Since the ratio Tc/ 11 is exponentially suppressed in weak coupling, one finds 1 3> e^i 3> at high baryon density. 
Note also the relation ~ lO^^ed/cA', the relative sizes of the critical regions are related to the types of color 
superconductor.^ Schematically shown in Fig. ^is a comparison of the sizes of the critical region of the scalar field 
and that of the gauge field for a type I superconductor. 

Let us briefly summarize the results obtained in this section. For a type I color superconductor, as realized in the 
high density region, gluon fluctuations give the dominant correction to the free energy. As we show in the next section, 
in the one-loop approximation these fluctuations change the order of the phase transition from second to first, and 
modify the critical temperature from Tc to T* . If the relative shift of the critical temperature (T* — Tc)/Tc is well 



^ This relation can be derived beyond the weak couphng approximation; the exact relation is = (l/256)e£j/^A- 
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outside the critical region dictated by Eqs. Hll() and (|12|l . the Gaussian approximation is consistent. This situation is 
quite analogous to the first order transition in type I metallic superconductor 

On the other hand, in a type II color superconductor, which may be realized in the low density region for Tc 
comparable to /i, scalar fluctuations are not at all negligible. Furthermore, the Gaussian approximation becomes 
highly questionable. A renormalization group analysis with an e expansion shows that, even without gauge fields, 
scalar fluctuations alone induce a first order transition in an SU{n) x SU{n) model with n > 3 Our model falls 
in this category when the coupling of the gauge field with the diquark condensate is neglected. A further complication 
for type II color superconductors is that the non-Abelian self-coupling of the gauge field may not be negligible, namely, 
the second condition Eq. (|13|l also becomes questionable at low densities. The non-Abelian coupling may change the 
order of the transition, as discussed in 0, but that problem is beyond the scope of this paper. 

The phase diagram in the (/9i,/92) plane implied by the GL free energy Eq. ^ is given in at the mean-field 
level without gauge fields. Once we take into account fiuctuations of the order-parameter and the gauge fields, we 
need to consider the phase structure in the four dimensional /?2, ht, as) space. Figure |2] shows its projection onto 
the two-dimensional cts) space with /3 = /3i = /32 = 3kt- The CFL phase in the weak coupling hmit lies in this 
reduced space (see the Appendix). 

The solid line in the figure is the boundary separating type I and type II color superconductivity, characterized by 
K = 1 in Eq. (|14|l . For the Gaussian approximation for the gauge fields to be reliable, 2as/7r should be smaller than 
1, Eq. (|13|l . Therefore, within the shaded region in the figure, the one- loop approximation taking into account only 
the gauge field is reliable for studying effects of thermal fluctuations. The dashed line in the figure shows the relation 
between and (3 (both functions of n) in weak coupling, as obtained from the formulas in the Appendix. The weak 
coupling regime is well within the shaded area. 




FIG. 2: The phase structure in the two-dimensional {f3,as) space. The solid line, characterized by k = 1, is the boundary 
between type I (right side) and type II (left side) superconducting behavior. The dashed line shows as as a function of /3, 
calculated in weak coupling. The loop expansion for gauge field fiuctuations is valid in the shaded area. 

In Fig. 131 we show k calculated in weak coupling, Eq. (|15|) . as a function of the baryon chemical potential. The 
dependences of jj, on Tc and as are taken from the weak coupling results in the Appendix. The figure indicates that 
K ^ 1 (type I superconductivity) is satisfied not only at high density but also at moderate densities, to the extent 
that one can rely on the extrapolation using the weak coupling formulas. 

III. FIRST ORDER TRANSITION INDUCED BY GAUGE FIELD 

In this section we assume a type I color superconductor and evaluate the free energy of the CFL and IS states up to 
one-loop order, taking into account gauge field, but not scalar field, fluctuations. The free energy difference between 



1 10 10^ 10^ 10^ 10^ 10*^ 
M- (GeV) 



FIG. 3: The parameter k calculated as a function of /i in the weak coupling approximation, Eq. (1151 . 

the superfluid and normal phases, _F'off(c?), in this approximation reads, in the Coulomb gauge, 

1 . . / 



FM = S{d) + - In det [^ Q^^Q^, ) - FMd - 0) 

= Sid) + ^ 2T^ / ^{In [k' + (m^)L] - In k'} 

= 5(d) [a2(mA)L - a3(mA)L +a4(nT-A)L] > (18) 



with 



We have included the two degrees of freedom associated with the polarization of the gauge fields in the Coulomb 
gauge. The high-momentum cutoff of the loop integral A, which sets the scale of the three-dimensional effective 
theory, is taken to be Tc, as discussed in the preceding section. The toa, the non-zero Meissner masses listed in Table 
^ are proportional to d, which here is a variational parameter determined by minimizing _Fcff in the above formula. 
In evaluating the integral we have used the expansion in terms of m/A = m/Tc. 

f(^) = / TTTTT log = / dAP 



|k|<A 



(2^)3 -V J Jo 7|k|<A(2^)3^^ + M2 

^m^-^m^ + -^m^ + 0{m'^). (20) 
27r2 67r 47r2A ^ ' ^ ' 



The expansion is well convergent for <^ Tc, corresponding to the condition Eq. (|12|l . However, this is not 
necessarily a good expansion for ruA ^ Tc, which is realized at the critical point of first order transition as we shall 
see in Sees. Ill A and III B. 

By setting m — nid in Eq. (|20(l . we can roughly estimate the energy contribution from diquark fluctuations. This 
equation also allows us to see the reason why k, introduced in the previous section, measures the relative importance 
of diquark and gauge field fluctuations in the free energy. At each order of the expansion of F^sid) in terms of d, 
the dominant contribution comes from the modes with largest mass. [Note that the energy scale entering Eq. H18|l . 
Tcwj^ A; can be regarded as the unit of thermal energy ~ T times the number density, ~ ^, of modes of mass 
md,A-] For the type I case, k — nid/rnA <C 1, which indicates that the gauge field fluctuations dominate over diquark 
field fluctuations. One can also see from Eq. 120() that the massless modes in Table I do not contribute to the free 
energy difference Feg. 

As we have discussed in the previous section and is easily seen by comparing Eqs. (|18|l and (Q, the integral of 
the gauge field fluctuation (A'^) with respect to leads precisely to the one- loop correction to the free energy. In 
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particular, the term proportional to 02 in Eq. (|18|) decreases the critical temperature of the second order transition 
since it is proportional to Cp' with a positive coefhcient. On the other hand, the term proportional to 03 drives the 
first order transition because it has a cubic structure with a negative coefficient. This is indeed the mechanism of 
the first order transition induced by thermal photons, first pointed out in 'll] in the context of the type I metallic 
superconductors. 

In dense QCD, Bailin and Love have previously discussed the first order transition induced by thermal giuons in 
two- flavor color superconductors 2] . Considering only the terms proportional to 02 and 03 , they concluded that the 
transition is strongly first order at high density. Furthermore, the strength of the first order transition grows as /i 
increases in their result [e.g., Eq. (4.100) in 0]. On the contrary, as we show below, the strong first order transition 
induced by the 03 term is substantially tamed by the consistent inclusion of the 04 term for both CFL and IS orderings. 
As a consequence, the first order transition becomes weaker with increasing density, in contrast to the conclusion of 



A. Free energy of the CFL state 



In this section, we show explicitly the one-loop free energy and the critical temperature of the first order transition 
T* for the CFL phase. If we use the CFL form in Eq. and the Meissner masses in Tabled the free energy Fcs{d) 
in Eq. (^HJ) reads 

(21) 



Feff(ci) = iad^ -qd^ +rd 



with 



3ad^ — iad^ + 8Ta2m\ ~ 



TTc 

3a -I — (32TrKTas 



qd^ 



8%/2r 



3/2 



and 



rd^ = 3(3/3i + (32)d'^ + 8Ta4m\ = 



3(3/3i+/32) 



128K^a? 



(22) 
(23) 

(24) 



Here, niA = {2KTy^^g\d\ is the Meissner mass in the CFL phase. 

As we have already mentioned, the effects of gauge fluctuations are threefold: First, they increase the size of the 
quadratic (d^) terms. This increase implies that thermal fluctuations tend to make the superconducting phase less 
energetically favorable. Second, we flnd a cubic term with a negative coefficient, —q, which favors the superconducting 
phase. Furthermore, the second-order transition found in mean-field theory turns into a first order one due to this term, 
independent of the magnitude of fj,. Finally, we find a positive correction to the quartic term from the fiuctuations. 
Like the correction to the quadratic term, this term acts against the superconducting phase. The sum of these three 
corrections leads to a first order transition which is significantly stronger than in a metallic superconductor, but is 
much weaker than that claimed in 2], where the quartic correction was neglected. 

For later convenience, we define a "renormalized" critical temperature T^. at which the quadratic term in Eq. H21|) 
vanishes, a = 0. (Note that the critical temperature at the tree level, Tc, has been defined through a — 0.) Then, 



= 1 



T2 32KTas 



T 



-T' 



T' 



37r 



aot'. 



(25) 
(26) 



The decrease of the renormalized critical temperature from the mean-field value {T'^ < Tc) is consistent with the fact 
that the contribution of the fluctuations to the quadratic term makes the superconducting state less favorable. 

On the other hand, the true critical temperature, T* , of the first order transition is the temperature where the 
free energy has two degenerate minima, d — and d = d* ^ 0. Since i^cff(O) = by definition, this implies that 
Fas{d*) = dF^s{d*)/dd* = is satisfied at T = T*. Thus we find 



T! 



= 1 



12aor ' 



d* 



2r' 



(27) 
(28) 
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The cubic term —q(P in the free energy, which tends to stabihze the superconducting phase, increases the critical 
temperature of the first order transition from T^, as seen in Eq. (|27|l . Thus the ratio T* /T'^ is a measure of the strength 
of first order transition induced by the thermal gluons. One may also define other measures such as the jump of the 
diquark condensate at T = T* relative to its T = value, d* /do, where do — d{T — 0). 

Using the weak coupling formulas for /3j and kt in the Appendix, we find the following estimates of the above 
quantities at high density: 



8l7r 



3anr 



g2 14C(3)a, ' 



T' 

— 

T, 



^ ~ 1 
T' 



5/2 



1/2 



d*_ 

do 



S^^3/2 / g xl/2 1 



10 \7C{3) 



1 

4^' 



H7 V 100 MeV 



1 GcV 



0.3k, 



(29) 
(30) 
(31) 
(32) 



where we have used the mean-field relation do = {n/e^)Tc with 7 being Euler's constant. 

In Fig.^Ja) we compare the ratio T^/Tc, T* /T^, and T* /Tc as functions of the baryon chemical potential fj,, with 
the weak coupling parameters of the Appendix. Note that T* is always smaller than Tc, 



T* 

— - 1 



as 



(33) 



In Fig. Elb), the jump of the order parameter at the critical point, d* /do, is shown as a function of /i. The jump is 
at most a few percent for /i > 1 GeV, but is much larger than that expected in type I metallic superconductors (see 
below). Also, one finds that the first order transition becomes weaker logarithmically as increases, and approaches 
a second order transition at /i = 00. As noted above, this is in contrast to the result of in which the first order 
transition becomes strong as increases. Such behavior is unreasonable because the coupling of gluons and diquarks 
becomes weak at high density due to the asymptotic freedom. As we have already discussed, the discrepancy between 
our result and that in originates from the fact that thermal corrections to the quartic term in the free energy were 
not taken into account in As we approach the baryon density close to that at the confinement-deconfinement 
transition, /i ~ 1 GeV, the fluctuations of the scalar field as well as the non-Abelian interactions of the gluons neglected 
in our calculation become important. Therefore the results shown in Figs. EIa,b) may be modified qualitatively in 
this region. Elucidating the super-to-normal transition in the low density region remains an interesting open question. 



1.02 



0.98 



0.95 





FIG. 4: (a) The ratio of three critical temperatures, Tl,/Tc (solid line), T* /T^ (dashed line) and T* /T'^ (dotted hne). (b) The 
jump of the diquark condensate at T = T* relative to its T = value; d* /do as functions of the baryon chemical potential fi 
in the CFL phase. Weak coupling values for /3i and kt are used. 
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It is instructive to compare the strength of the present first-order transition with that of type I metallic superconduc- 
tors. In the latter case, one finds an extremely weak first-order transition: [T* — T^)/T^ ~ (Ep /Tc)"^ (Ep /m)'^^'^ 
10~^ Ep is the electron Fermi energy {Ep = p|./2to), m is the electron mass and a = /A-k is the electromag- 

netic fine structure constant. Major differences from the case of massless color superconductivity are the presence of 
a small factor Ep/m and the small coupling constant. Note that in the metallic case, the thermal photon correction 
to the quartic term in the GL free energy divided by the corresponding mean-field value is 0[a'^{Ep/TcY{Ep/m)], 
which is the same order with or even smaller than unity. In contrast, thermal gluon fluctuations in the weak-coupling 
color superconductors dominate the quartic term [see Eq. (|24ll ]. This is the reason why in the metallic case the shift 
(T* — T^)/T^ involves such a high power, a'^, compared with only for color superconductivity. 

Let us now discuss the reliability of the first order phase transition obtained here in weak coupling, from the point 
of view of the critical region examined in Sec. IIIBI The first inequalities in Eqs. Hll|l and \i'2.l can be interpreted 
as conditions for the size of T — Tc, or alternatively as conditions for the size of d; if T* — Tc or d* is very small, 
critical fluctuations are not negligible and one cannot trust the result of the one-loop approximation. In the case of 
weak coupling, t* = \T* — Tc\/Tc ~ lO^^a^ from Eq. (|33|l : in addition TcZ/i <C 1. Therefore, the conditions given in 
Eqs. ((TTjl and lfTB|l . t* ^ e-A,d, are well satisfied, and the critical temperature of the first order transition is outside the 
critical region of the diquark and gauge fluctuations. Also, by substituting d* in Eq. H31|l into the first inequalities in 
Eqs. (|ll|l and H12|) . one finds the conditions as/ fi ^ 1.6 tt^ and ^ ""^/B, which should be satisfied for the effective 
three-dimensional approach to be valid. They are, in fact, well satisfied, insofar as the couplings stay in the shaded 
region in Fig. [2| 

B. Free energy of the IS state 

The analysis of the IS state is similar to that of the CFL state. The free energy in the one-loop approximation 
becomes 



Feff(d) = ad^ - qd^ + rd'^, 



(34) 



with 



ad^ 



qd^ 



rd" 



ad -f Ta; 



T 



Here, uia and 



iPi + I32)d'^ + Tai (4m^ 
are the Meissner masses. Table 



2 

3 + 

4 

- TO ; 



TTc f 327TKxa 



208 
~9~ 



2 2 
Krra„ 



2 72 

KTg d 



(4 components), 
(1 component). 



(35) 
(36) 
(37) 

(38) 
(39) 

/3i < /32 < 



If the global minimum of the mean- field theory is in the IS state, that is, the parameters /3i, /?2 satisfy 
(as shown in Fig. 1 of 0), the effective potential in the mean field approximation yields a second order transition to 
the normal phase, and in the one- loop approximation, a first order transition. The renormalized temperature, T^, the 
critical temperature of the first order transition, T*, and order parameter, d* , at the minimum of the free energy are 
calculated as before: 



T' 

^ C 

Tc 
d* 



T2 32KTa, 
ao 37r 



C 

T' 



Aa^r ' 



2r" 



(40) 
(41) 



The relation, T*/Tc < 1, is also satisfied in this case. Since some of to^'s in Table U are negative in the weak 
coupling (/3i = /32 > 0), IS ordering is unstable against scalar fluctuations and decays into CFL ordering at high 
densities. This result is consistent with the result in Ref. Q where the comparison of the free energy between the 
IS and CFL orderings is made. As found from such comparison, the IS state is not in a local minimum when the 
CFL state is in the energy minimum. Nevertheless, the above formulas are relevant at finite temperature where the 
unlocking transition from CFL to IS ordering takes place due to the effect of the strange-quark mass, TOs jHES- 
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IV. SUMMARY 

In this paper we have studied the effect of thermal fluctuations of diquarks and the gluons on the superconducting- 
to-normal state phase transition in a three-flavor color superconductor. For this purpose, we adopted the Ginzburg- 
Landau free energy in three-spatial dimensions. Although the phase transition of this model is second order in the 
mean-field approximation for coupling constants near those in weak coupling, it can be turned into a first order 
transition either by the thermal fluctuations of the scalar diquark field, or the giuon gauge field near the critical point. 

The relative importance of these two types of fluctuations is controlled by k, the ratio of the masses of the scalar 
field and the gluon just below the critical temperature; k is also the Ginzburg-Landau parameter that differentiates 
type I and type II color superconductors. In the high density regime where the weak coupling approximation is valid, 
we find that the system is type I and gauge fluctuations dominate over scalar fluctuations. 

After evaluating the size of the critical region, outside of which the one-loop approximation is a reasonable ap- 
proximation, we calculate the one-loop correction to the free energy from the thermal gluons in a type I color 
superconductor. The transition to color superconductivity becomes first order due to the induced cubic term in the 
GL free energy, which is similar to the case of the type I metallic superconductors pTl |. 

The strength of the first phase transition can be characterized by quantities such as the change of the critical 
temperature from its tree-level value, and the jump of the order parameter at the critical point. They indicate that 
the first order transition weakens with increasing baryon density. This behavior, which is quite reasonable in the 
sense that gluonic corrections are suppressed by as in weak coupling, is in sharp contrast to that found for two-flavor 
color superconductivity in 0]. The difference stems from the fact that one needs to take into account not only the 
cubic term of the order parameter (which strengthens the first order transition) but also thermal corrections to the 
quartic term (which suppresses the first order transition) to obtain a consistent result. Since the Ginzburg-Landau free 
energy has an intrinsic spatial cutoff scale A~^, which is of order the size of the diquark (~ T~^), such a correction 
to the quartic term is an inevitable consequence of the three-dimensional effective theory. We also find that the 
critical temperature of the first order transition, T* , is always lower than the Tc of the second order transition in the 
mean-field approximation. 

Our general considerations in this paper for a type I superconductor are valid insofar as the parameters in the 
Ginzburg-Landau free energy stay in the shaded area in Fig. [3 which corresponds to the high density region. On the 
other hand, in the low density, strong coupling region, not only scalar fluctuations but also non-Abelian interactions 
among thermal gluons are not negligible. Anti-quark pairing and the non-perturbative running of as at low momentum 
are also not negligible at low density 23]. These effects may change the nature of the phase transition at low density 
from that realized at high densities. Lattice simulations of the SU{3)c x SU{3)l+r x U{1)b sigma model -I- SU{3)c 
gauge field introduced in the present paper would be a good starting point to analyze the phase structure in the 
strong coupling region. Furthermore, the strange quark mass plays an important role in the unlocking transition from 
the CFL state to the IS state. 
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APPENDIX: ASYMPTOTIC VALUES OF THE COUPLING CONSTANTS 

In this paper, we use Tc calculated in finite temperature perturbation theory in the normal phase [2^: 



TT 



67r 

= 4^ " (33 - 2Nf) ln[A./(3AQCD)] ' ^ ' 

where b'^/tt 0.57 and Aqcd for A^/ = 3 is taken to be 200 MeV. 

The parameters in Eq. as a function of and /i (baryon chemical potential) at asymptotically high densities, 
which were calculated from the finite temperature weak coupling gap equation in Refs. ,^,SJ, are 

a = AN{^i/i)\n(^\ (44) 



7C(3) 
8(7rT,)2 

1 //iN2 



/?i=/32 = 3KT = ^7f^A^(/i/3), (45) 



with the zeta function C(3) = 1.2020 • • • . 
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